Eigenstate thermalization within isolated spin-chain systems 
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The thermalization phenomenon and many-body quantum statistical properties are studied on the example of 
several observables in isolated spin-chain systems, both integrable and generic non-integrable ones. While diag- 
onal matrix elements for non-integrable models comply with the eigenstate thermalization hypothesis (ETH), the 
integrable systems show evident deviations and similarity to properties of noninteracting many-fermion models. 
The finite-size scaling reveals that the crossover between two regimes is given by a scale closely related to the 
scattering length. Low-frequency off-diagonal matrix elements related to d.c. transport quantities in a generic 
system also follow the behavior analogous to the ETH, however unrelated to the one of diagonal elements. 
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Many-body quantum (MBQ) systems and models have 
been extensively studied in the last decades in connection with 
novel materials, offering a fresh view on the fundamentals 
and the interpretation of statistical mechanics. The systematic 
analysis of the phenomena of thermalization and the limita- 
tions of a statistical treatment within isolated MBQ systems 
have been recently motivated by experiments on cold atoms 
in optical lattices, revealing very slow relaxation to thermal 
equilibrium [1, 2], but as well by prototype integrable MBQ 
systems as the one-dimensional Heisenberg model realized in 
real materials [3]. 

Specific for lattice MBQ systems discussed in the above 
connection is (in contrast to single-body quantum systems) 
the exponential growth of the Hilbert space and the number of 
eigenstates N st with the lattice size L. Here, one of the fun- 
damental questions is to what extent even a single eigenstate 
or a single chosen initial wave-function could be the repre- 
sentative of the canonical ensemble average within the given 
system, both for static and dynamical quantities. For generic 
MBQ systems one of the central statements is the eigenstate 
thermalization hypothesis (ETH) [4, 5] that for a few-body 
observable A diagonal matrix elements (ME) A aa at a given 
energy show only exponentially (in L) small deviations from 
the average, being a smooth function of the energy only. Since 
at the same time the off-diagonal ME are as well exponentially 
small, the time-average of the observable is determined by di- 
agonal terms only. Therefore for any initial wave-function 
|^o) with a small energy uncertainty the long-time average is 
also equal to the thermal average, this being the general condi- 
tion for the quantum thermalization process [6]. We note that 
such a hypothesis is also underlying some numerical methods 
for the calculation of finite-temperature properties, in partic- 
ular the microcanonical Lanczos method (MCLM) [7, 8] for 
T > static and dynamical properties of MBQ lattice sys- 
tems. It seems also evident that the ETH is intimately re- 
lated to general properties of eigenenergy spectra, i.e. level 
statistics and dynamics in generic MBQ systems, which reveal 
Wigner-Dyson level statistics with the origin in level repulsion 
and analogy to random matrix spectra [9, 10]. 

The deviations from the ETH and normal thermalization 



have been detected is several directions. The hypothesis is not 
obeyed in integrable MBQ systems [6, 1 1-14], although some 
observables can still thermalize, i.e., approach the equilib- 
rium (canonical ensemble average) value, in particular if the 
Gibbs statistical ensemble is generalized to include all local 
conserved quantities in this case [11, 14]. The thermalization 
can become very slow and the validity of ETH can become re- 
stricted if an initial state is far from equilibrium [12, 15-17] as 
relevant for sudden quenches in cold-atom systems. The latter 
question is intimately related to the deviation from integrabil- 
ity [13] and the size of isolated MBQ systems [6, 15, 17]. On 
the other hand, the ETH does not resolve the question of the 
relation to off-diagonal ME (even in generic non-integrable 
systems) which are, e.g., relevant for transport properties and 
dissipation in the d.c. limit [18-20]. 

In this Letter we study the validity of the ETH and ther- 
malization within a quantum spin-chain system in one dimen- 
sion (ID), i.e., the antiferromagnetic and anisotropic S =1/2 
Heisenberg model (AHM), including integrable (I) and non- 
integrable (NI) cases. While we confirm in the generic NI case 
the ETH for diagonal ME of several local observables, we find 
large deviations and fluctuations for the I case. In particular, 
we show that the spread of diagonal ME can be qualitatively 
and even quantitatively understood from the XX model equiv- 
alent to the model of noninteracting fermions. With the aim 
to resolve the problem of the breakdown of the ETH in fi- 
nite systems we perform the finite-size scaling in NI systems 
revealing that the crossover from the I regime to the ETH- 
consistent behavior is determined by a single scale L*, coin- 
ciding with the transport scattering length. Another finding is 
that the low-frequency off-diagonal and diagonal ME are not 
universally related even in NI systems, hence the ETH does 
not directly address the low-w dynamics and the d.c. transport 
quantities, and the generalization of the ETH is necessary. 

As the prototype model we study in the following the 
anisotropic 5 = 1/2 Heisenberg model on a chain with L 
sites and periodic boundary conditions, 
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where S" (a — x, y, z) are spin 5 = 1/2 operators at site i 
and A represents the anisotropy. The nearest-neighbor XXZ 
model is an integrable one and we introduce the next-nearest- 
neighbor zz-interaction with A 2 7^ in order to break its in- 
tegrability. It should be reminded that the Hamiltonian (1) 
can be mapped on a t-V-W model of interacting spinless 
fermions. A consequence of the integrability at A2 = is the 
existence of a macroscopic number of conserved local quanti- 
ties and operators Q n ,n= 1 — L commuting with the Hamil- 
tonian, [Q n > H] = 0. For the ID AHM a nontrivial example is 
Q3 = J E representing the energy current and leading directly 
to its non-decaying behavior [21, 22] and dissipationless ther- 
mal conductivity [3]. 

In order to study eigenstate and ME properties of the AHM 
we choose some simple local (but q = 0) quantities involving 
only few sites. Evident candidates are nontrivial quantities 
involving n = 2 sites, where we consider the "kinetic" energy 
H km containing the first two terms in Eq. (1) and the spin 
current J 3 , while for a representative of n = 3 operators we 
consider the energy-current J E (not including the A2 term), 
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Figure 1 . (Color online) Distribution of diagonal matrix elements of 
(a) spin current J s , (b) energy current J E , and (c) kinetic energy 
H km vs. energy E for the integrable model A2 = (l.h.s.) and the 
non-integrable model A2 = 0.5 (r.h.s.). In all cases: A = 0.5, 
L = 20, M = -1, and k = 2tv/L. 
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The choice is motivated by different properties of the consid- 
ered operators. While J E is a strictly conserved quantity for 
the I case, J s is not, but still leads to dissipationless (non- 
decaying) spin transport. Both are current operators with ME 
distributed around the ensemble average (J^m) = 0. On the 
other hand, H^ n has not such a specific property. In the fol- 
lowing we present results reachable via the exact diagonal- 
ization of the model, Eq. (1), on chains up to L = 20. The 
total spin = M is fixed to M = —1 (in order to avoid 
"particle-hole" symmetry) while we consider both, the repre- 
sentative sector with wavevector k — 2ir/L and the whole 
fc-average as well. 

First, we present results for the distribution of diagonal ME, 
i.e., J* n , J E n , and H^", as they arise varying eigenenergies 
E = E n . In Fig. 1 we show corresponding 2D plots obtained 
within the gapless regime (A = 0.5) and for the magneti- 
zation M = — 1 (due to "particle-hole" symmetry J 3 n van- 
ishes at M = 0). Figure 1 reveals an evident difference be- 
tween the NI example with A2 = 0.5 and the I case with 
A2 = 0. All quantities show for the NI example a narrow 
distribution around the average A(E) = (A nn )(E) with the 
width (a£(E)) 2 = (A* n ) - (A nn ) 2 exponentially dependent 
on the system size L [6]. 

On contrary, for the I case distributions are much wider with 
a weaker size dependence, clearly not obeying the ETH. The 
distribution for J s and J E is intimately related to the anoma- 
lous T > spin and energy-current stiffness (without degen- 



eracies), respectively, for the I model [18, 19, 22], 

D s,E (T) = tf_J2e-^\J^\ 2 , (3) 

n 

where (3 s = /3j E = f3 2 with (3 = 1/T. It is evident that 
the existence of D S ' E (T > 0) > implies that currents as 
J S ' E do not thermalize to their thermal average (J S ' E ) = 0. 
In particular, their correlation functions do not decay to zero, 
(J s ' E (t — > oo)J s ' E ) 7^ 0, and their time evolution depends 
crucially on the ensemble of initial states. The same appears 
to be the case for H km , although a physical interpretation is 
less familiar. With values of D S ' E (T) known from the Bethe 
Ansatz [21], and moreover for the energy-current stiffness 
D E (T — > 00) obtained easily via the high-T expansion, one 
can evaluate the distribution widths <j^' E (E) oc \f~L. 

Since analogous quantities to stiffness are not known in 
general, one can use in the gapless regime (A < 1) as a semi- 
quantitative guide results for the A = XX model. The lat- 
ter can be mapped to the model of non-interacting fermions, 
H = J^k e k n k, £fc = Jcosfc, being trivially integrable with 
all rife = 0, 1 as constants of motion, with corresponding cur- 
rents J S ' E = J2k(^> £ k)vknk- The calculation of a^' E (E) 
at fixed magnetization M = J2k( n k — 1/2) averaged over 
energies E is for L — > 00 equivalent to the grand-canonical 
averaging in the limit /3 — > yielding for the unpolarized case 
N = L/2: 0% = Jy/Z/V8 and a E = J 2 VI/\/32. On the 
other hand, instead of H kln (being within the XY limit equal 
to H) one can treat in an analogous way the complementary 
potential term H A with the result erf = JAvT/4 [23]. We 
note that the above estimates for the widths cr^ represent well 
the numerical results in Fig. 1 for the I case with A > 0. 

Next we investigate the crossover from an I to a NI sys- 
tem obeying the ETH. In a finite system fluctuations a% = 



Figure 2. (Color online) Finite-size dependence of the diagonal ME 
fluctuations of J\ J E , and -ff kin , respectively, for A = 0.5, M = 
— 1 and different A2. Fluctuations are evaluated within E = [—1,1], 
In (b) the exact sum rule is indicated (solid line). Inset in (a): Curve 
for A2 = 0.5, illustrating the exponential decrease with L. 

<7j / VX with a = (s, E, kin) are expected to decrease by in- 
troducing the NI perturbation A2 7^ 0. In Fig. 2 we present 
corresponding results obtained for different A2 = — 0.5 and 
sizes L = 8 — 20. In order to reduce the influence of the 
energy window, we evaluate the fluctuations cr^ in the range 
E = [—1,1] and average over all fc-sectors. For the I case 
A2 = the 1/L-scaling indicates finite values a$(L —t 00). 
This coincides with the well defined and nontrivial D S ' E (T — > 
00). In particular D E (T — > 00)/ f3 s ' E and af can be related 
to the high-T sum rule (erf ) 2 = (1 + 2A 2 )/32 [21]. This is, 
however, not the case for the NI case A2 ^ 0. Here, there 
is an evident decrease with L and crossover to an exponen- 
tial decrease with L, i.e., ETH-consistent behavior above the 
crossover scale L > L*. L* crucially depends on the per- 
turbation strength A2 but apparently is quite universal for all 
quantities, at least for those considered here. 

In the case of currents the "thermalization length" L* may 
be plausibly interpreted in terms of the transport mean free 
path. The latter can be determined by a standard hydrody- 
namic relation, l/(q 2 T>) 3> 1/7 [24], involving the diffusion 
constant V and the current scattering rate 7. Identifying the 
mean free path as L* = ir/q then yields 

L* = n^V/^. (4) 

In the case of the spin current, using for A = A2 = 0.5 
the known quantitative values V s — 2.1 and 7 s = 0.23 at 



Figure 3. (Color online) Probability distribution for the real part of 
the off-diagonal ME of J s , J E , and ff Mn , respectively, for a NI 
model with A = A 2 = 0.5, L = 18, and M = -1, For com- 
parison, Gaussian functions are indicated (dashed curves). 

/3 [20], one finds L* ~ 10. This value turns out to agree 
well with the scale observed in the inset of Fig. 2. Moreover, 
7 s — > as A2 — > is consistent with a diverging L*. 

Finally, let us address the relation between off-diagonal and 
diagonal ME. Since for the I system the behavior can be very 
singular [20], we concentrate on the generic NI cases satisfy- 
ing the ETH. In Fig. 3 we present the probability distribution 
of off-diagonal ME Re and Reff^, evaluated for A = 
0.5 and A 2 = 0.5 in the energy window E = [-SE/2, SE/2] 
with various SE. Resulting distributions do clearly not de- 
pend on SE and appear to be Gaussian with the width <r™ d , 
again in analogy to ETH exponentially decreasing with L, It 
is a nontrivial observation that both off-diagonal and diagonal 
ME follow the same scaling with L, as shown in the inset of 
Fig. 4. It is therefore important and well defined to investigate 
the ratio of off-diagonal and diagonal ME fluctuations 

r a (E)=a« l (E)/a2(E). (5) 

Results for the spin and energy current are presented in Fig. 4, 
shown vs. E for A2 = 0.5 and A = 0.5,1.0. They indi- 
cate that r a (E) is not universal (depends on a and model pa- 
rameters) and smoothly varies with E, but most important is 
the independence of L. We can conclude that for the cases 
considered here r a are quite far from the predictions of the 
random-matrix theory [9, 19] implying generally r = 2. On 
the other hand, the ratio still remains within an order of mag- 
nitude in contrast to the I case where in the gapless regime the 
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Figure 4. (Color online) Ratio between off-diagonal and diagonal 
ME fluctuations for J s and J E , respectively, for a NI model with 
A = 0.5, 1, A 2 = 0.5, L = 20, M = -1, and fc-average. Inset in 
(a): Finite-size dependence of r a (E = 0). 

ratio appears to vanish leaving finite only diagonal ME [20]. 

The above observation becomes relevant in the evaluation 
of d.c. transport quantities, which are within linear response 
theory related to the low-w absorption [25], e.g., the spin con- 
ductivity (diffusivity) and thermal conductivity, respectively, 
are in analogy to Eq. (3), 

C Q H = ^ ]T e-^"|J«J 2 <^ - E m + E n ) , (6) 

where the d.c. limit should be considered as C§ = C a (u> — > 
0) and can be expressed as 

C °=^rJ z- pE P 2 {E)K A {E)] 2 dE, (7) 

where p(E) is the MBQ density of states. From our analysis it 
follows that in general a^ A (E) cannot be represented by diag- 
onal crj (E), although the qualitative behavior appears closely 
related. Note that for the case of J s diagonal ME can be also 
expressed as the sensitivity of many-body levels to a fictitious 
flux <fr (or boundary conditions), i.e. J* n oc dE n /d(j>, and the 
latter relation has been previously employed to evaluate the 
d.c. transport in, e.g., disordered systems [10, 26]. 

Let us in conclusion summarize our results, which appear 
to be generic beyond spin-chain systems. The behavior of the 
considered NI systems we find consistent with the ETH for 
all considered quantities. If we consider the time evolution of 
an observable, it can be in terms of (finite-system) eigenstates 
represented as 

A(t) = (tt(t)|A|*(t)> = ]T \c n \ 2 A nn + 

n 

+ J2 <c m e^"- B ™)% mi . (8) 

n^rrt 



In a system obeying ETH, the off-diagonal contribution van- 
ishes for long times t — >• oo, due to the exponential smallness 
of off-diagonal ME as well due to dephasing [6]. If the initial 
state | "Jo) is a microcanonical one with a narrow distribution 
SE [with (SE) 2 = J2 n \c n \ 2 {E n - E) 2 ], and due to ETH 
A nn ~ (A)(E), the first term leads to the microcanonical av- 
erage A(t) = (A)(E) in a large system coinciding with the 
canonical thermodynamic al average at a finite T > 0, where 
E(T) = E. Such a scenario is then consistent with the "nor- 
mal" quantum thermalization. 

In an I spin chain the distribution of diagonal ME is large, 
the long-time average [still neglecting off-diagonal terms in 
Eq. (8)] in general depends on |^o) an d corresponding c„, 
even for a small energy uncertainty SE. In order to satisfy 
A(t — > oo) = (A) one needs assumptions on the distribution 
of coefficients c„. E.g., in a large enough system randomly 
chosen c n would plausibly be adequate. In fact, the numerical 
MCLM method for the evaluation of T > properties [7, 8], 
based on the microcanonical states and the Lanczos proce- 
dure, contains such a choice achieved by random sampling. 
Hence, a random microcanonical state in a large MBQ sys- 
tem would mostly obey the thermalization process. Still, this 
is not at all the case for particular states as, e.g., reached by 
(strong) quenching in an I system, but as well not in a generic 
system [13, 17] since the initial state after the quench is not 
necessarily the microcanonical one with small SE. 

Analyzing the extent of the validity of the ETH and ther- 
malization in a finite-size MBQ system, we find effectively 
that perturbed I systems beyond the crossover length L* be- 
have as generic NI ones. Since in a "normal" spin system 
only total spin and energy are conserved, one can design two 
relevant diffusion scales and plausibly the largest would de- 
termine L*, which then appears to dominate the scaling of all 
quantities, as shown in Fig. 3. The understanding and the de- 
termination of L* is evidently an important theoretical goal, 
relevant also for experiments dealing with systems close to 
integrability [3, 27]. 

The ETH addresses thermalization and statistical descrip- 
tion of static quantities in MBQ systems, with the behavior 
determined by diagonal ME. On the other hand, d.c. transport 
quantities and low-w dynamics involve only off-diagonal ME. 
We note that in a generic system, properties analogous to the 
ETH can be defined for off-diagonal ME close in energy, in 
particular obeying the Gaussian distribution and exponential 
dependence on size. Also, the relation between diagonal and 
off-diagonal ME is independent of size L, but still the ratio is 
not universal. In this sense, our results show that for such con- 
siderations the generalization of the ETH is needed but also is 
straightforward, and it can include the response to weak ex- 
ternal fields and dissipation phenomena in MBQ systems. 
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